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We provide an update of the determination of the electromagnetic self-energy contribution to 
M p — M„ based upon Cottingham's Formula. A technical oversight in the literature is uncovered: 
the application of the Cottingham Formula requires the use of a subtracted dispersion integral; 
an argument to evade the subtraction function was presented; the argument was based on false 
£■ — assumptions about the scaling violations of the parton model, a point first mentioned by J. C. 

Collins. We elucidate this point and utilize low-energy effective theory to relate the unknown 
subtraction function to the nucleon isovector magnetic polarizability. This allows us to provide 



(SI the first reliable determination of 8M? = 1.30(03)(47) MeV [1]. 



X 



The 30 International Symposium on Lattice Field Theory - Lattice 2012, 
June 24-29, 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http : //pos. sissa. it/ 



Cottingham Formula 



Andre Walker-Loud 



1. Introduction 

Given only electrostatic forces, a natural prediction is M p > M„ yet in nature [2] 

M p -M n = -1.29333217(42) MeV. (1.1) 

Before we knew of QCD, there were many attempts to reconcile this apparent discrepancy, see 
Ref. [3] for a review. We now know the Standard Model has two sources of isospin breaking: the 
electromagnetic couplings of the light quarks and their mass parameters. The contributions to the 
nucleon mass splitting from these two sources are roughly equal in magnitude but opposite in sign. 
While the net result is well known, Eq. (1.1), our ability to disentangle the contribution from these 
two sources remains relatively poorly constrained. 

An understanding of Eq. (1.1) from first principles is desired. For example, a determination of 
the electromagnetic contribution to the nucleon mass splitting will allow for an independent means 
of determining mj — m u [4]. At leading order (LO) in isospin breaking, we can write 

8M p - n = M p - M„ = 8M r + 8M m "- m " . (1.2) 

The quark mass operator is needed to renormalize electromagnetic self-energy contribution, en- 
tangling these effects beyond LO. In general, electromagnetic effects can not be unambiguously 
separated from hadronic effects so one must formulate a prescription to disentangle them [5]. At 
LO, lattice QCD can be used to reliably compute the contribution from m t \ — m u . A weighted, 
uncorrected average from 3 independent lattice calculations [6] yields, 1 

Sm lqcd" = -2.53(40) MeV . (1.3) 

Calculating 8M y with lattice QCD and QED is much more challenging, due to the disparate length 
scales relevant for the two theories. There is an alternative method for computing 8M y known 
as the Cottingham formula [8], which relates the electromagnetic self-energy to an integral of the 
forward Compton scattering tensor. In 1975, Gasser and Leutwyler provided the determination [9] 

8M r = 0.76(30) MeV . (1.4) 

If we instead subtract the lattice QCD determination from experiment, we arrive at 

8M p - n - 8M'l d QC ™" = 1.24(40) MeV . (1.5) 

While the central values in these two determinations are quite different, they are only one-sigma 
discrepant. A more precise determination of either 8M y or 8M md ~ m « is clearly desirable. We 
report on an effort to improve the electromagnetic self-energy contribution using the Cottingham 
Formula [1]. 



2. The Cottingham Formula 

The Cottingham Formula intuitively relates the electromagnetic self-energy to forward Comp- 
ton Scattering (see Fig. 1) 



'There is an additional determination in Ref. [7] which we have not yet included. 
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(a) (b) 

Figure 1: The Cottingham Formula intuitively relates the forward Compton Scattering (a) to the ff(af. s .) 
electromagnetic self-energy (b). 

The subscript R on the integral reminds us that the self-energy contains a logarithmic divergence 
and must be renormalized [10]. The spin-averaged forward Compton Scattering tensor is given by 

Tnv(p,q) = ^Lf d A xe^ x {pa\T{j^x)J v (0)}\pa) 

-jr-gnv)Ti(q - q )+[Pn-qn-^rJ [Pv-qv-fJ — ^ — < 2 - 2 ) 

Cottingham showed that by performing the Wick rotation q Q — > iv and then a variable transforma- 
tion Q 2 = q 2 + v 2 , the nucleon self-energy can be related to the experimentally measured structure 
functions; 

where 2 

t£ = -3T1 (iv.fi 2 ) + ( 1 ~y) Ti(iv,Q 2 ). (2.4) 

One uses fixed-fi 2 dispersion integrals to determine the scalar functions 7}(iv, Q 2 ) in terms of their 
experimentally measured absorptive (imaginary) parts. It is known that T 2 satisfies an unsubtracted 
dispersion integral while T\ requires one subtraction [11]. These scalar functions are crossing 
symmetric 7}(— v, Q 2 ) = 7}(v, Q 2 ), thus given by 

v 2 r°° 2v' 

7\(v,!2 2 ) = ?i (0, Q 2 ) + - JdV (v02((v/)2 _ y2) Zfari (V + ie, Q 2 ) (2.5) 
1 f°° 2v' 

T 2 (v,Q 2 ) = =- dv' 2lmT 2 (v' + ie,Q 2 ) (2.6) 

in Jv, h [y ) - v 

The +ie in the argument indicates the function is evaluated just above the cut on the positive real 
axis (see Fig. 2). The absorptive parts are given in terms of the well known nucleon structure 
functions 

2Imr 1 (v,e 2 ) = 27r^i(v,e 2 ), 2lmT 2 (v,Q 2 )=2x^F 2 (v,Q 2 ), (2.7) 

to which we have implicitly included the isolated elastic nucleon pole. 2 Except in the low and high 
Q 2 limits, the subtraction function Ti(0,Q 2 ) can not be simply related to measured cross sections, 
complicating the determination of 8M r . 

2 The elastic pole is isolated because we are working to leading order in QED, so there is no Bremsstrahlung radia- 
tion. The inelastic cut begins at the pion production threshold. 
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Figure 2: Analytic structure of scalar functions 7}(v,g 2 ) 
contour used to evaluate them. 



in the complex v plane for fixed Q , and the 



The presence of the unknown subtraction function severely impacts the ability to determine 
8M y precisely. In Ref. [9], an argument to evade the subtraction function based on the parton 
model was presented. However, as was first noted in Ref. [10], the argument was based on false 
assumptions about the scaling violations of the Callan-Gross relation [12]. To understand the issue, 
we first recall that the decomposition of the Compton tensor, Eq. (2.2) is not unique. Another 
common decomposition is 



Tnv(p,q) = (guv- ( * } ^J q 2 ti(q°,-q 2 



p-q i . \ Ap-q) 2 
PuPv — 2" KP^qv + Pvq^i ) h 2~ Suv 



q 2 t 2 (q ,-q 2 ) 
M 2 



(2.8) 



What are the advantages of this choice over Eq. (2.2)? For a point particle, such as the electron, at 
leading order in QED, t\ = 0. In general 

2lmt l (v,Q 2 ) = 2n ^ V [2xF l (x,Q 2 )-F 2 (x,Q 2 )} where x = . (2.9) 

In the deep-inelastic (DIS) limit 2xF\(x) —F2(x) = 0. In Ref. [9], it was assumed the corrections to 
this Callan-Gross relation were power law, 

HUx) 

IxF^-Fj^t-^ 1 , (2.10) 

as predicted by the parton model. This extra suppression in v would allow for an evasion of 
the unknown subtraction function. However, it is known in QCD that the corrections are only 
suppressed by the strong coupling constant [13]; 

2*F 1 (*,<2 2 ) -F 2 (x,Q 2 ) = - 3 ^°^1f 2 (x,Q 2 ) (2.11) 

In the Regge limit (fixed Q 2 , v — > °°,x — > 0), it is known emperically [1 1] 

limF'^" oc => limlmff~" oc <fi (2.12) 

x^0 *->0 
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such that in the unsubtracted dispersion integral for ff - ", the contour at infinity does not vanish but 
in fact diverges, see Fig. 2. One is lead to conclude that the unknown subtraction function can not 
be evaded and the evaluation of 8M y in Ref. [9] is not correct. 

Using either representation for T^ v , Eqs. (2.2) or (2.8), performing a once-subtracted disper- 
sion integral at v = 0, one arrives at the same answer 

8M y = 8M eI + 8M inel + 8M sub + 8M ct . (2.13) 

The renormalization is detailed in Ref. [10]. The operator product expansion (OPE) was used to 
connect the ultraviolet behavior of the integrand of Eq. (2.3) with operators in the QCD+QED 
Lagrangian. After canceling the log-divergence, there is a residual finite contribution, 8M ct , which 
can be estimated using arguments similar to Naive Dimensional Analysis [ 14] ; 

tar * - ( i) m ~ 5S iflf \"l CM* 

An VAg/ 9m (p\uu + dd\p) 

where o k n = m(p\iiu + dd\p) /2M ~ 45 MeVandthe quark masses are given by 2m = m u + m c i and 
28 = m d - m u . Taking A 2 , = 2 GeV 2 and A 2 = 100 GeV 2 , one can estimate \8M ct \ < 0.02 MeV. 

The most difficult contribution to determine precisely is that from the unknown subtraction 
function 3 

•sub _ JCC / j/->2t i(\ rt 2\ 



= - I- d ® Tl (°> Q ) • ^ 2 - 15 ^ 



o 



The OPE constrains the high Q 2 behavior, XirciQi^ T\ (0, Q 2 ) °< l/Q 2 and low-energy effective 
theory constrains the low Q 2 behavior [15] 

Mm li (0, Q 2 ) = 2k(2 + k) - Q 2 1 \ [(1 + k) 2 v 2 m ~r 2 E ]+^~ 2M^ } + 0(Q A ) , (2-16) 

where K is the anomalous magnetic moment, r £ ( M j is the electric (magnetic) charge-radius and 
j8m is the magnetic polarizbility. To connect the high and low Q 2 behavior, one must resort to 
modeling. All but the last term of Eq. (2.16) are recognized as the leading terms of elastic form 
factors, suggesting a resummation. While one cannot derive the resummed formula from first 
principles, it can be computed from the Born graphs of nucleon-Compton scattering with the full 
elastic form factors inserted at the vertices [16]. The low Q 2 limit of Ti(0,Q 2 ) is now known 
to €?(<2 4 ) from heavy baryon #PT and the expansion of the resummed Born contributions agrees 
precisely through this order [17]. Using the resummed elastic terms, one arrives at 4 

Tr\0,Q 2 )^2G 2 M {Q 2 )-2FKQ 2 ) + 2 -^^Q 2 [l-^^ , (2.17) 

with M 2 = -ft/ 3M(AnF z m % f F 92 MeV and the isoscalar magnetic mom ent u s ~ 0.88. A re- 

cent review of low-energy Compton scattering yields the value " = — 1.0± 1.0 x 10 -4 fm 3 [19]. 
If the two inelastic terms were of opposite sign, we could re-sum them to a dipole form factor as in 



3 The same unknown subtraction function contributes to (muonic) hydrogen energy levels [15, 16, 17, 18]. 
4 We thank J. McGovern for providing the isovector formula to 0(Q 4 ) from Ref. [17]. 
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Ref. [17]. However, the large uncertainty in even the sign of j3^~" makes this impractical to imple- 
ment in a numerical evaluation. To be conservative, we can simply multiply the leading inelastic 
contribution by a standard dipole form factor such that the overall contribution obeys the high Q 2 
scaling. This leaves us with a natural separation of the subtraction term 



SM: 



sub . 



-3a 



dQ 2 [G 2 M (Q 2 )-F 2 (Q 2 )], 8MZ b , 



-3fljT" 

87T 



dQ 2 Q 2 



llln 



ml + Q 2 



(2.18) 



3. Evaluation of 8M^ 



Details of the numerical evaluation can be found in Ref. [1]. The elastic and inelastic contri- 
butions are given by the expressions 



8M e 



a 



^IdQ 2 f3^G 2 M + [G 2 E -2z el G 2 M ] 



1 + T d 



(l + T^ 2 -^ 2 - 



(3.1) 



8A? 



inel 



a 



^dQ 2 



4xM Jo Q 



dW 2 



3Fi(v,6 2 ) 
M 

, F 2 (v,Q 



t 3 / 2 -t^/TTt + ^/t/2 



+ 



T 

3/2 



r 3/2. 



(3.2) 



where x d = Q 2 /4M 2 , % = V 2 /Q 2 , W 2 = M 2 + 2Mv - Q 2 and W 2 h = (M + m n ) 2 . The value of 
Aq = 2 GeV 2 was used. The elastic terms are evaluated with the Kelly parameterization of the 
elastic form factors [20] and the inelastic contributions are evaluated from parameterizations of 
the resonance [21] and scaling regions [22]. The subtraction terms are estimated from Eq. (2.18), 
conservatively taking rnk = 0.71 GeV 2 . For the isovector self-energy, we find 



8M el = 1 .39(02) MeV , SM'f 

SM? ub 



8M" iel =0.057(16) MeV : 



inel 



-0.62(02) MeV, 
+0.47(47) MeV, 



(3.3) 



for the total 



8M y = 1.30(03) (47) MeV. 



(3.4) 



4. Concluding Remarks 

A precise determination of the electromagnetic self-energy contribution to M p — M n remains 
an outstanding theoretical challenge. We have updated the Cottingham determination of 8M y using 
modern theoretical knowledge and experimental determinations of the nucleon structure functions, 
Eq. (3.4). In the process, a technical oversight in the old determination was uncovered invalidating 
the result. The limiting factor in a precise determination comes from modeling the behavior of the 
unknown subtraction function between the low and high Q 2 limits where it is known. Presently, the 
large uncertainty can be traced directly to the poorly determined isovector magnetic polarizability. 
Once this quantity is precisely known, an improved understanding of Tf~ n (0, Q 2 ) will be required 
to make progress. 
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